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ABSTRACT 

We present results from an optical spectroscopic investigation of the massive 
binary system [L72] LH 54-425 in the LH 54 OB association in the Large Mag- 
ellanic Cloud. We revise the ephemeris of [L72] LH 54-425 and find an orbital 
period of 2.24741 ± 0.00004 days. We find spectral types of 03 V for the primary 
and 05 V for the secondary. Wc made a combined solution of the radial velocities 
and previously published V-band photometry to determine the inclination for two 
system configurations, i = 52I3 degrees for the configuration of the secondary 
star being more tidally distorted and z = 55 ± 1 degrees for the primary as the 
more tidally distorted star. We argue that the latter case is more probable, and 
this solution yields masses and radii of Mi = 47± 2 Mq and Ri = 11.4 ± 0.1 Rq 
for the primary, and M2 = 28 ± 1 Mq and i?2 = 8.1 ± 0.1 Rq for the secondary. 
Our analysis places LH 54-425 amongst the most massive stars known. Based on 
the position of the two stars plotted on a theoretical HR diagram, we find the 
age of the system to be ~1.5 Myr. 

Subject headings: binaries: spectroscopic - stars: early-type - stars: fundamental 
parameters - stars: individual ([L72] LH 54-425) 



1. Introduction 

One of the most elusive of important parameters in stellar astrophysics is the mass of 
stars. Stars in binary or multiple systems, and specifically in double-lined spectroscopic 



^ Based on observations with the CTIO SMARTS Consortium l.S-meter telescope and the 2.5-meter 
DuPont telescope at Las Campanas. 
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binaries, offer us a cliance to explore tfieir masses witfi one important caveat: we must know 
ttie orbital inclination of the system. The inclination is the key that allows the resolution of 
the spectroscopic results into individual masses for the members of the system. Individual 
stellar masses are obtained via the serendipitous orientation of the system for eclipses or, 
at the very least, ellipsoidal light variations. Particularly rare are 0-type main-sequence 
systems, owing partly to their short lifespans and partly to the steepness of the initial mass 
function. Reported here is the first of several systems to be investigated with the goal of 
determining the masses, radii, and distances to a number of 0-star binaries. Mass data of this 
kind are needed to help resolve the mass discre pancy problem between the evolutionary a nd 
spectroscopic mass estimates for massive stars (IRepolust et al.l 12004 iMassey et al.ll2005l ). 



The e arliest type star in the Large Magella nic Cloud LH 54 OB association is [L72] 



L H 54-425 flHill et ai. 



by lOey fc Smedleyl ( 



19941) . also known as L54S-4 (IOeylll996l ). Its spectral typ e was classified 



19981 ) as 04 III(f*). Its binary nature was unknown until Ostrov ( 20021 ) 



noted the ellipsoidal variations in the V^-band photometry and obtained six spectroscopic 
observations of the system. The light and radial velocity curves from that study were used 
to estimate the masses of the two stars at 100 Mq and 50 Mq. The lack of phase coverage 
for the radial velocity data did not allow for a more rigorous analysis of physical parameters. 

We present here an analysis of the radial velocities from optical spectra obtained between 
2003 and 2006 (§2). In §3 we discuss how we extracted radial velocities from the observations 
and solved for the orbital parameters. We describe tomographic reconstructions of the 
individual component spectra in §4. The next section (§5) covers our combined light curve 
and radial velocity solution, as well as our determination of the inclination of the system 
for two system configurations. We conclude with a discussion of the system's evolutionary 
status, fundamental parameters, and a comparison of the two model fits to the system in §6. 



2. Observations 

The first sets of spectroscopic data were obtained from 2003 through 2006 with the 
Cerro Tololo Inter- American Observatory (CTIO) 1.5-m telescope and RC spectrograph 
operated by the SMARTS Consortium. We used grating #47 (831 lines per mm, 8000 A blaze 
wavelength in the first order Littrow configuration) in second order (SMARTS configuration 
47/IIb) together with a BG39 or CUSO4 order blocking filter. The detector was a Loral 
1200 X 800 CCD. This arrangement produced spectra covering the range 4058 - 4732 A with 
a resolving power of i? = A/AA ^ 3550. Exposure times were 1800 seconds, in order to reach 
a moderate signal-to-noise ratio (S/N ^ 70). Each observation was bracketed with HeAr 
comparison spectra for wavelength calibration, and numerous bias and flat field spectra were 
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obtained each night. We obtained a total of 48 spectra of the target with this configuration. 

In addition to the CTIO data, four observations were obtained with the echelle spec- 
trograph on the 2.5-m DuPont telescope at Las Campanas Observatory binned 2x2 with a 
resolving power of ~ 27000 (from measurements of arc spectra). Exposures were combined 
where there were two, and the spectral flux was normalized. Orders were merged and cosmic 
rays removed. For continuity of analysis, these data were clipped and re-sampled to match 
the wavelength and lower resolution CTIO data. 

The resulting spectra were extracted and wavelength calibrated using standard routines 
m IRAF@. All spectra were then rectified to a unit continuum via fitting line-free regions and 
transformed to a common heliocentric wavelength grid in log A increments. 



Radial Velocities and Orbital Elements 



Radial velocities were measured using a template-fitting scheme (iGies et al.ll2002l ) for 
the five lines B.6 A4101, He II A4200, H7 A4340, He II A4542, and He II A4686. These were the 
only prominent lines in the spectrum of each star, and no single line showed any systematic 
shift from the others. The spectra show no evidence for wind effects in the lines used to derive 
the orbital solution. No clear emission features and no asymmetries of absorption features 
were seen in the hydrogen lines or t he He II A4686 line. The basic procedure followed to 
obtain radial velocities is outlined in iBoyajian et al.l (120071 ). 



The spectrum from HJD 2,452,718.5370 shows well separated spectral features from 
each star in the He II AA4542, 4686 lines. These lines were used to derive matching template 
spectra for the primary and secondary components. Templat e spectra were obtained from 
the grid of 0-type stellar models from lLanz fc Hubenyl (120031 ). which are based on the line 
blanketed, non-LTE, plane-p arallel, hydros t atic atmosphere code T LUSTY and the radi ative 
transfer code SYNSPEC dHubenvlEoSsI : iHubenv fc Lan3 Il995l : iHubenv et aDll998h . In 
order to find matching templates, approximations were initially used for the temperatures, 
gravities, projected rotational velocities, and flux contributions from each star. These initial 
parameters for model templates were then checked after velocity analysis by studying the 
properties of the tomographically reconstructed spectra of the components (§4). 



Further inputs for the template fitting scheme are preliminary estimates of the velocities 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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for the primary and secondary. The He II A4542 hne for each star was fit using Gaussian 
functions with the splot routine and deblend option in IRAF. Radial velocity shifts were 
then calculated for all spectra with well separated lines, and these were used to obtain initial 
orbital parameters. The preliminary velocity from the initial orbital solution gave a starting 
point to perform a non-linear, least squares fit of the composite profiles with the template 
spectra and calculate the shifts for each star. The five values for radial velocity from each 
set of five lines in an individual spectrum were then averaged to obtain a mean value, and 
the standard deviation of the mean value was calculated. Each of these values are listed in 
Tabled] for the primary and secondary stars, as well as the orbital phase for each observation 
and the observed minus calculated values for each point. Zero phase is defined by the time 
of inferior conjunction of the primary star, Tic,i. 

We made an orbital fit of the radia l veloc ity data using the non-linear, least-squares 



fitting program of iMorbey fc Brosterhud (119741 ). Equal weights were assigned to each data 



point, because all the spectra were comparable in S/N ratio, spectral coverage, resolving 
power, and quality of radial velocities derived. It should be noted that for massive binaries, 
discrepancies in systemic velocities between components may be attributed to differences in 
their expanding atmospheres and/or differences in the shapes of template spectral lines, so 
we fit systemic velocities for each component. Fits with non-zero eccentricities were explored, 
none of which gave as good a fit as a circular orbital solution. Fits were also made for subsets 
of measurements near conjunction (phases near 0.0 and 0.5) and quadrature (phases near 
0.25 and 0.75). These tests revealed that the important spectroscopic elements, such as the 
velocity semi- amplitudes, did not vary by more than 3 a from the best fits for any subset. 
Our results for the orbit using spectroscopic data only are listed in column 3 of Table [2l The 
radial velocity measurements and final orbital velocity curves (§5) are shown in Figure [H 



4. Tomographic Reconstruction 



We used a Doppler tomography algorithm (iBagnuolo et al.lll994l ) to separate the pri- 
mary and secondary spectra of [L72] LH 54-425. This is an iterative method that uses the 
48 observed composite spectra from CTIO, their velocity shifts, and an assumed monochro- 
matic fiux ratio {F2/F1) to derive the individual component spectra. We explored a range 
of fiux ratios to arrive at a value that best matched the line depths in the reconstructions 
with those in model spectra. Figure [2] shows the reconstructed spectra for the primary and 
secondary, as well as synthetic spectra just below each. Few lines are present in our spectra. 
However, the relative strength of the He I A4471 in the secondary is the major difference 
between the two spectra and indicates the cooler temperature of the secondary. 
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These reconstructed spectra were fit with T LUSTY/ SY N SPEC model synthesis spec- 
tra (see §3). These matches allow for estimates of stellar parameters listed in Table [31 The 
width of the He II A4542 line was used to estimate the projected rotational velocity V sin i 
for each star by comparing the reconstructed and model profiles for a grid of test V sini 
values. Matches were made to the H7 A4340 and He I A4471 lines to obtain the temper- 
ature and surface gravity estimates. Finally, we used the ratios of the reconstructed to 
synthetic spectral line depths to estimate that the monochromatic flux ratio {F2/F1) in our 
blue spectra is 0.45 ± 0.10. 

A visual inspe c tion o f our reconstructed spectra match best the 03 V ex ample in Figure 



6 oflWalborn et al.l (|2002| ) for the primary and the 05 V example in Figure 7 of lWalborn et al. 



(I2OO2I ) , and these are the spectral types we list in Table [31 Unfo rtunately, our wavelen gth 
coverage does not include the N IV A4058 feature suggested in IWalborn et al.l (|2002| ) for 
spectral classification of the hottest 0-stars. The N III A4634-42 feature is another good 
diagnostic for classification of early 0-type stars, but is too weak in our spectra for r neasure- 
ment. We can use the equivalent width ratio of He I A4471/He II A4542 calibrated in iMathys 
(119881) for each of our reconstructed spectra to obtain a spectral class. This measurement 
results in a spectral class for the primary of 03 and of 05.5 for the secondary. We used our 
best matched effective temperatures from t he T LUSTY/ SYNSP EC models combined with 
the Teff versus spectral type calibration of [Martins et al.l (|2005[ ) to estimate spectral types 
of 03 V for the primary and 05 V for the secondary. These sp ectral types and effec tive 
temperatures are also consistent with those shown in Figure 1 of [Mokiem et al.[ (120071 ) for 
LMC 0-type dwarfs. 



5. Combined Radial Velocity and Light Curve Solution 



Light curve data were taken from the \^-band observations listed in Table 1 of [Ostroy 
([20021). New orbital parara eters were then found using the Eclipsing Light Curve (ELC) code 



(lOrosz fc Hauschildtl [20001 ) . which fits the radial velocity and light curves simultaneously. 
The resultant light curve is shown in Figure [3l and the radi al veloc i ty cur ves are shown in 
Figure m along with the radial velocity measurements from [Ostrovl (|2002| ). which were not 
used for any fit in this paper. ELC's genetic optimizer was used to explore the parameter 
space and was given wide ranges for each value of period, epoch of inferior conjunction of 
the primary Tic,i, inclination, mass ratio, primary velocity semi-amplitude, and Roche- lobe 
filling factor for each star. The Roche-lobe filling factor is defined as the ratio of the radius 
of the star toward the in ner Lagrangian Point (Li) to the distance to Li from the center of 
the star, / = Xpoint/a^Li (lOrosz fc Hauschildtl [2OOOI ): our fits show that our stars do not fill 
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their Roche lobe (see Table H]). A by-product of the exploration of parameter space by ELC 
was a determination of systemic velocities for each component. 

The temperatures of each star were set and not fit because these were well constrained 
by analysis of the reconstructed spectra discussed in §4. As with the fits of only radial 
velocities discussed in §3, circular orbital solutions proved best. 

To estimate the 1, 2, and 3 a uncertainties as sociated with the fitted parameters, we fol- 



lowed the method discussed in lOrosz et al.l (120021 ). The values of the seven fitted parameters 



(plus the two systemic velocities) listed above were varied in the calculation of ~3.6xl0^ 
light and radial velocity curves. This procedure resulted in thorough sampling near the Xmm 
point, making the multi-dimensional surface well populated. This allowed for projection of 
the surface as a function of each parameter or fitted value of interest. The lowest 
is then found for each parameter of interest, and the 1, 2, and 3 a uncertainties may be 
estimated by the regions where x^ < Xmm +1; +4; and +9, respectively. The uncertainties 
listed in the tables are 1 a uncertainties. 

Special attent i on wa s paid to the inclination of the system. The well-sampled light 



curve from lOstrovl (120021 ) allows for an exploration of the important parameters that cause 
the ellipsoidal variations in this system: inclination and Roche-lobe filling factors for each 
star. Ellipsoidal light variations occur in systems where one or both of the components are 
tidally distorted due to the proximity of its companion. At different orbital phases, the 
observer sees varying amounts of light from the system, based upon the interplay between 
the amount of tidal distortion, the observed cross-section of the two stars, and the inclination 
of the system. As was expected, we found that for lower inclinations, i.e., viewing the system 
more face-on, the sizes and tidal distortions of the stars needed to become larger in order 
to match the modulation of the light curve. Thus, restricting the stars to radii within their 
Roche lobes ultimately limits the lowest acceptable inclination, while the lack of observed 
eclipses establishes the maximum inclination. 

Again, u sing the genet ic algorithm in ELC, we explored the range of inclination values 



suggested in lOstrovl (120021 ) and allowed the fill factors for each star to vary over a wide 
range while keeping the other parameters fixed. The results of this analysis are plotted as 
contours of 1, 2, and 3 a confidence intervals in Figures H] and |5l In the best overall fit, 
the estimated filling factors are larger for the secondary, so most of the ellipsoidal variation 
is due to the tidal distortion of the secondary (Fig. H]). It is interesting to note that the 
contours suggest a rather steep surface toward higher inclinations, indicating a hard upper 
limit for the inclination of 56° at the 3 a level. Lower inclinations are not as probable, and 
fill factors for the secondary get too large to match the light curve as the secondary star 
approaches a fill factor of 1.0, thus completely filling its Roche-lobe. These characteristics 
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help us to constrain the inchnation of the system to be z = 52_3 degrees. The corresponding 
masses and other parameters are hsted in Table HI Included are values for Vgync sin i, the 
projected rotational velocity assuming synchronous rotation, i?efr, the effective radius for a 
sphere of the same volume, and the sizes of the stars along different directions. 

The orbital parameters for the best fit are listed in Table [2l The spectroscopic solution 
and the combined solution from ELC match very well. It should be noted that the parameters 
in Table H] with asymmetric uncertainties are a result of the asymmetry in the uncertainties 
for the inclination. 

We also explored the possibility of the primary being the more tidally distorted star. To 
do this we ran the simulations again, keeping the same parameters from Table [2] and varying 
only the filling factors and inclination. Thus, the radial velocity curve does not change and 
the light curve is fit to roughly the same degree of accuracy as before, so the fits in Figures 
[U and [3] need not be changed. The contour plot for this scenario, again with confidence 
intervals of 1, 2, and 3 a, is shown in Figure [51 The fit gives an inclination of i = 55 ± 1 
degrees. This corresponds to slightly different values of masses and radii, listed alongside 
the overall best fit in Table HI The small range in inclination leads to a small range in the 
masses and radii of the components of the system. This restricted fit has a reduced of 
2.69 versus 2.54 for the case where the secondary is more tidally distorted. 

A contact binary is another possible configuration for the system. We did not explore 
this possibility because in a contact system, the sum of the values of Vsini and the sum 
of the velocity semi-amplitudes would be comparable. For our fits of [L72] LH 54-425 (see 
Tables [2] and [3]) , the sums of these values differ by at least 300 km s~^. 

6. Discussion and Conclusions 

Which of our two models is most consistent with the available data? The answer to this 
question lies in comparing the consequences of each solution. 

To begin the comparison, we used the temperatures and radii output from ELC to 
obtain an estimate of the monochromatic flux ratio F2/F1. The case where the secondary 
is more tidally distorted yields a model flux ratio of 0.70 tooi, while the primary as the 
more tidally distorted star gives a flux ratio of 0.46 ± 0.02, in agreement with the value of 
0.45 ± 0.05 derived from the tomographic reconstructions. 



Theory suggests (lEkstrom et al.ll2008l ) that for massive stars, the main sequence radius 
scales as i? ~ M°-^. Applying this as another consistency check, we can substitute velocity 
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semi-amplitudes to estimate the expected ratio of radii, R2/R1 ~ {K1/K2) which, from 
our values in Table [2] is 0.707 ± 0.10. The ratio of effective radii estimated from the output 
of ELC (Table H]) is ~ 0.881o;q5 for the fit with the secondary as the more tidally distorted 
and ~ 0.71 ±0.010 for the primary being more tidally distorted. Once again the fit with the 
primary as the more tidally distorted matches theory more closely. 



Ignoring the tidal distortions, we can estimate the luminosity by applying the Stefan- 
Boltzmann equation L = AnR'^aT^g, and combining this with our temperature measure- 
ments, place each component in a t heoretical H-R diagra. m. This is presented in Figure El 
along with evolutionary tracks from ISchaerer et al.l (Il993l ) for stars of varying mass with a 
metallicity of Z = 0.008 which is appropriate for the LMC The 1 a uncertainty regions are 
consistent with an age of ~1.5 Myr for the model consisting of the primary as the more tidally 
distorted star, and both stars appear to be co-evolutionary. The tidally distorted secondary 
fit appears a bit older, roughly matching the 2 Myr isochrone. This fit also indicates that 
the secondary is more evolved than the primary and is overluminous for its derived mass. 
Therefore, it seems the fit with the primary being more tidally distorted is more consistent 
with evolutionary track s. These approx i mate ages are also consistent within uncertainties 

of ~2-3 Myr for the LH 54 OB association. 
, even the oldest age of ~2 Myr puts [L72] 



with the estimate from lOev fc Smedlev 
According to Table 1 of 



Schaerer et al. 



(1199^ 



LH 54-425 in the part of its life when it is still burning hydrogen, and it will continue to do 
so for another ~ 1.5-2 Myr. 

For the next comparison, we used the luminosities calculated above to obtain a bolo met- 
ric ab solute magnitude, assuming a solar bolometric absolute magnitude of 4.74 mag (I Cox 
2000). Nex t, we applied a bolom etric correction for each star based on our derived spectral 
types from lMartins et al.l (120051 ) to obtain absolute V magnitudes. We then combined these 
magnitudes to obtain an absolute V magnitude for the system. Finally, we arrived at a 
distance modulus for the system by using th e maximum magnitude of the system from the 
photometric data and E{B — V) = 0.10 mag (lOey &: Smedleylll998l ). The resulting distance 
modulus from the fit with the more tidally distorted secondary is 18.65 mag while the dis- 
tance modulus for the tidally distorted primary fit is 18.55 mag. Both of these numbers are 
consistent with the HS T Key Project distance modulus to the LMC of 18.50 ± 0.10 mag 
( iFreedman et al.ll200ll ). with the tidally distorted primary fit being slightly better. Note 
that each of the contour plots in Figures H] and [5] seem to have a best fit valley in the 
as a function of inclination and filling factor. As one travels along this valley to higher 
inclinations, the dimensions of the system get smaller, including the sizes of the two stars. 
This makes the stars intrinsically fainter, and they must lie at a closer distance to match the 
observed V photometry. In contrast, the lower inclinations mean the stars become bigger 
and lie at greater distances. The overall agreement in distance modulus between the model 
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results and those for the LMC indicate that the final fill out estimates are reliable. 



The two models are most different in the V^ync sin i values output from ELC. The 
observed V sin i values from the reconstructed spectra are 197 ± 5 km s^^ for the primary 
and 182 ± 8 km s~^ for the secondary. The ELC fit for the more distorted secondary gives 
values of 197 ± 5 km s~^ and 173 ± 7 km s~^ for the primary and secondary, respectively, 
while the ELC fit for the more distorted primary gives 209 ± 5 km s~^ and 148 ± 6 km s~^. 
Both cases agree within uncertainties for the primary's line broadening, but the predicted 
Vsini for the secondary is too small for the case where the primary is the more tidally 
distorted star. However, due to the youth of the system, our assumption of synchronous 
rotation may not be appropriate. There may not have yet been enough time for tidal forces 
to synchronize the spin with the orbit. More ELC fits were performed assuming a faster 
rotation for the secondary and these show that the same physical parameters are obtained 
for a non- synchronous rotation scenario. 

Given the above tests, we think the fit with the primary as the more tidally distorted 
star is the more likely configuration. For this case, the masses and radii are Mi = 47±2 Mq 
and Ri = 11.4 ± 0.1 Rq for the primary and M2 = 28 ± 1 Mq and R2 = 8.1 ± 0.1 Rq for 
sec ondary. These masses are significantly less than the 100 Mq and 50 Mq values obtained 
by lOstrovl (120021 ). Part of the reason lies in the data set we obtained for radi al veloc i ties o f 
the system. As is shown in Figure [T|, two secondary velocity values derived by lOstrovl (120021 ) 
near phase 0.25 are the ones that drove the estimate for the velocity semi-amplitude up, 
and conseque ntly the mass ratio down. This, in addition to the lower inclination adopted by 
Ostrovl (I2OO2I ). leads to the differences in ma sses between o ur two analyses. Our temperatures 



are also slightly lower than those given in lOstrovl (120021 . see hi s Table 3), due t o his use 



of a different, older temperature versus spectral- type cali bration (lAUer et al. 



compare our primary star wit h the primary in R136-38 (IMassey et al.ll2002l ). which is also 



1982h . We 



can 



an 03 V. iMassey et al.l (120021 ) derive a mass for their 03 V star of 56.9 ± 0.6 Mq, in closer 
agreement with our value of 47 Mq than that of 100 Mq found by lOstrovl (120021 ) . Our derived 
radius for the primary in LH 5 4-425 is slightly larger than that for the primary in R136-38 
of 9.3 Rq (IMassey et al.ll2002l ). This may be due to the extrem e youth of R136-38 , which 
appears to be near the zero-age main sequence (see Fig. 6 in IMassey et al.ll2002l ). while 
the position of [L72] LH 54-425 in the H-R diagram (our Fig. [HD indicates the primary is 
slightly evolve d. These numbers are consistent with the fact that the LH 54 OB association, 
at ~2-3 M vr (lOev fc Smedleylll998r) . is slightly older than the R136 cluster with an age of 
~l-2 Myr jMassev fc Huntei]ll998h . 
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Table 1. LH 54-425 Radial Velocity Measurements 



Date Orbital Vi ai (O - C)i V2 (Ji {O - C)i 

(HJD-2,400,000) Phase (km s"!) (km s"!) (km s"!) (km s"!) (km s"!) (km s"!) 



52714.5496 


0.080 


393.0 


6.7 


-4.1 


126.4 


20.7 


-16.2 


52714.5756 


0.092 


398.5 


12.9 


-11.1 


72.8 


49.6 


-47.5 


52714.6009 


0.103 


414.3 


4.9 


-6.7 


90.9 


42.1 


-9.2 


52714.6312 


0.116 


430.9 


11.7 


-3.1 


91.2 


14.6 


14.0 


52714.6586 


0.129 


431.4 


15.6 


-13.7 


39.1 


38.6 


-18.4 


52715.5251 


0.514 


269.6 


9.7 


-11.8 


405.9 


16.7 


57.1 


52715.5487 


0.525 


263.0 


9.1 


-4.8 


404.0 


23.1 


31.4 


52715.5725 


0.535 


231.3 


5.7 


-23.8 


423.6 


11.5 


28.5 


52715.5971 


0.546 


242.3 


11.1 


0.9 


384.8 


9.1 


-34.4 


52715.6216 


0.557 


213.4 


5.0 


-14.5 


418.4 


14.3 


-24.4 


52716.5492 


0.970 


245.8 


19.2 


-16.0 


389.8 


20.5 


4.4 


52716.5722 


0.980 


234.2 


17.4 


-40.4 


398.1 


17.8 


35.6 


52716.5952 


0.990 


252.9 


24.1 


-34.6 


354.5 


27.1 


15.3 


52716.6182 


0.000 


343.9 


31.6 


43.4 


231.1 


35.6 


-84.7 


52717.5139 


0.399 


443.8 


14.0 


23.5 


149.6 


39.9 


43.9 


52717.5370 


0.409 


419.3 


16.6 


10.1 


147.8 


39.3 


22.5 


52717.5642 


0.421 


404.9 


10.2 


8.8 


171.2 


21.5 


22.8 


52717.5878 


0.432 


387.5 


15.0 


3.5 


183.0 


23.3 


13.3 


52717.6122 


0.443 


378.2 


16.3 


6.8 


198.6 


17.0 


6.8 


52718.5139 


0.844 


139.1 


6.4 


5.6 


621.2 


55.7 


7.8 


52718.5370 


0.854 


137.6 


20.6 


-3.3 


595.3 


28.4 


-5.0 


52718.5600 


0.864 


142.0 


9.3 


-7.0 


563.3 


26.1 


-22.8 


52718.5831 


0.875 


151.1 


10.2 


-6.6 


575.4 


26.5 


4.6 


52718.6060 


0.885 


157.9 


7.4 


-9.2 


536.5 


18.1 


-17.8 


53017.5789 


0.915 


152.5 


12.4 


-44.9 


469.5 


23.3 


-31.0 


53017.7698 


0.000 


323.6 


25.6 


23.8 


259.5 


36.7 


-57.5 


53018.5366 


0.341 


499.9 


22.4 


30.8 


67.6 


46.7 


48.4 


53019.5645 


0.798 


102.2 


11.9 


-6.7 


654.1 


12.7 


-1.9 


53019.7591 


0.785 


160.2 


21.6 


-6.8 


555.3 


18.5 


0.8 


53020.5382 


0.232 


515.8 


11.4 


17.4 


-10.0 


42.2 


24.9 


53021.5305 


0.673 


149.7 


15.4 


28.0 


662.8 


8.9 


32.1 


53021.7435 


0.768 


104.3 


6.8 


3.2 


670.6 


3.7 


1.5 


53022.5223 


0.114 


437.1 


23.5 


5.0 


135.2 


47.4 


54.7 


53023.5216 


0.559 


249.2 


10.4 


23.7 


419.2 


8.9 


-27.8 


53027.5315 


0.344 


441.2 


14.8 


-26.2 


25.6 


8.6 


3.5 


53028.5271 


0.786 


118.3 


17.8 


13.2 


674.8 


17.5 


12.3 


53028.7270 


0.875 


165.0 


8.0 


6.9 


591.8 


37.4 


21.7 


53029.5269 


0.231 


518.6 


9.8 


20.2 


19.6 


32.4 


54.4 


53030.5260 


0.676 


131.5 


10.3 


11.4 


631.7 


15.5 


-1.9 


53030.7169 


0.761 


96.7 


10.5 


-3.6 


683.5 


20.2 


13.2 


53031.5200 


0.118 


475.7 


22.6 


40.1 


118.6 


20.8 


44.2 


53032.5246 


0.565 


234.0 


18.9 


15.5 


487.5 


48.1 


28.3 


53032.7225 


0.653 


146.3 


3.9 


12.1 


569.6 


21.2 


-38.7 


53033.5391 


0.017 


346.1 


6.3 


25.4 


209.5 


20.0 


-69.8 


53710.7543^ 


0.347 


539.9 


39.3 


75.4 


31.7 


33.4 


5.4 
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Table 1 — Continued 



Date Orbital Vi 


cri (O - C)i 


V2 0-2 (O - C)2 


(HJD-2,400,000) Phase (km s"! 


) (km s-l) (km s-i) 


(km s^^) (km s"-') (km s"-*-) 


53711.7167=- 0.776 77.1 


52.7 -29.5 


688.2 79.4 23..6 


53748.5493 0.165 458.6 


11.7 -13.4 


14.7 21.9 4.2 


53752.6857 0.006 304.4 


25.3 -2.3 


282.6 31.5 -22.0 


53753.5655 0.397 454.1 


15.9 32.0 


73.8 23.6 -28.8 


53755.7542 0.371 438.2 


9.6 -8.5 


-3.4 20.9 -62.5 


53765.6205'' 0.761 55.1 


21.2 -49.6 


623.4 64.7 -44.4 


53766.6131=- 0.203 526.3 


43.6 40.6 


-54.9 78.2 -37.6 


"Data from the echelle spectrograph on 


the 2.5 m DuPont telescope at Las Campanas. 


Table 2. Circular Orbital Solutions for [L72] LH 54-425 


Element 


Combined Solution 


Spectroscopic Solution 


P (days) 


2.24741 ± 0.00004 


2.24746 ± 0.00010 


TiCi (HJD-2,400,000) . . . 


53029.007 ± 0.003 


53029.016 ± 0.007 


Ki (km s-^) 


201.6 ± 3.8 


210.8 ± 3.0 


K2 (km s-i) 


359.1 ± 4.5 


350.8 ± 5.2 


7i (km s-^) 


299.7 ± 0.6 


303.1 ± 2.0 


72 (km s~^) 


316.7 ± 2.2 


317.7 ± 3.3 


rms (Primary) (km s~^) . 


19.4 


21.9 


rms (Secondary) (km s~^) 


33.2 


32.8 


rms (Photometry) (mag) . 


0.007 





Table 3. Tomographic Spectral Reconstruction Parameters for [L72] LH 54-425 



Parameter 


Primary 


Secondary 


Spectral Type 


03 V 


05 V 


Teff (kK) 


45 ± 1 


41 ± 1 


log g (cgs) 


4.0 ± 0.2 


4.0 ± 0.2 


V sin i (km s~"^) 


197 ± 5 


182 ± 8 


F2/F1 (blue) 


0.45 ± 0.05 
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Table 4. ELC Model Parameters for [L72] LH 54-425 



Parameter 


Distorted Secondary 


Distorted Primary 




Primary 


Secondary 


Primary 


Secondary 


Inclination (deg) 


52 


+2 
-3 


55 ± 1 


M(Mo) 


53 tl 


32 


47 ± 2 


28 ± 1 


ReS {Rq) 


11.0 


q 7 +1-0 
-0.2 


11.4 ± 0.1 


8.1 ± 0.1 


Rpole^ {Rq) 


10.7 t'oi 


q 2 +0-3 
^•^ -0.2 


10.9 ± 0.1 


7.8 ± 0.1 


-fi^point*^ (-^0) 


1 1 7 +0.4 
-0.2 


11 1 +°-^ 
^^■^ -0.2 


12.3 ± 0.2 


8.6 ± 0.1 


V^ync sin i (km s""*^) 


197 ± 5 


173 ± 7 


209 ± 5 


148 ± 6 


log g (cgs) 


4.08tHl 


3.941H? 


4.00 ± 0.02 


4.07 ± 0.01 


Filling Factor 


U.UU _Q Q2 


0.80 


r, 79 +0.03 
' ^ -0.02 


64 

•J-"^ -0.02 


fltot (Rq) 


31.6 


+1.0 
-0.6 


30.4 


± 0.4 


F^/F^ (blue) 


'-'•"-'-0.01 


0.46 ± 0.02 



^Polar radius. 

''Radius toward the inner Lagrangian point. 
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Orbital Phase 



Fig. 1. — Radial velocity curves for [L72] LH 54-425. Primary radial velocities are shown by 
filled dots and secondary radial velocities are represented by filled triangles with associated 
uncertainties shown as line segments for both. The open circles and triangles represent the 
echelle data from L as Canipanas. The open diamonds and squares are the radial velocity 
measurements from lOstrovl (120021 ) for the primary and secondary, respectively. The Ostrov 
data were not used in model fits. The solid line is the best combined fit solution for the 
primary, and the dotted line is the same for the secondary. 
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Fig. 2. — Tomographic reconstructions of the components of [L72] LH 54-425 based on 48 
spectra obtained from 2003 through 2006 at CTIO. Plotted from top to bottom are hne 
i dentifications with ver tical tick marks, the primary spectrum, the model primary spectrum 
(ILanz fc Hubenyll2003l ). the secondary spectrum, and the model secondary spectrum. The 
stellar parameters for the model spectra are given in Table [31 
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Fig. 3. — \^-band light curve for [L72] LH 54-425. These data were taken from lOstrovl (120021 ) 
and are presented here in phase according to our best combined solution. The model is the 
solid line and the data are represented by filled dots with V uncertainties shown by line 
segments. Phase zero corresponds to inferior co ni unction of the primary star (which differs 
by 0.5 phase from that adopted by IOstrovll2002l ). 



-19- 




Fig. 4. — surface contours of the residuals from the combined velocity and light curve 
data as a function of orbital inchnation and secondary Roche-lobe fiUing factor of the best 
fit solution. The best fit position is represented by the "x" inside the l-cr contour. The range 
of filling factor is quite large, while the inclination is reasonably well constrained. 
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Fig. 5. — surface contours of the combined fit residuals as a function of the orbital 
inclination and Roche-lobe filling factor of the primary star for the case where the primary 
is the more tidally distorted star. The best fit position is represented by the "x" inside the 
1-a contour. The parameter space is more confined than for the best fit, and the contours 
are smaller around the best inclination for this configuration. 
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Fig. 6. — An H-R diagram showing the location of the primary star {filled circle) and 
secondary star {filled triangle) of [L72] LH 54-425 for our fit where the primary is the 
more tidally distorted star. The open symbols represent the fit where the secondary is the 
more tidally distorted star . Also plotted are evolutionary tracks for stars of various masses 
from ISchaerer et al.l (119931 ) for an LMC metallicity. The boxes around the filled data points 
correspond to the 1-cr uncertainties in the derived values of T^s and R from the fits where 
the primary sta r is the more tidally distorted star. The vertical dashed lines correspond to 
isochrones from lLejeune fc Schaererl (120011 ). also for an LMC metallicity, of 1, ~2 and ~3.2 
Myr going from left to right. The positions of the two stars are consistent with an age of 
~1.5 Myr. 



